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Abstract 

QCD predicts soft radiation patterns that are particularly simple for VF + jet pro- 
duction. We demonstrate how these patterns can be used to distinguish between the 
parton- level subprocesses probabilistically on an event-by-event basis. As a test of our 
method we demonstrate correlations between the soft radiation and the radiation inside 
the outgoing jet. 

1. Introduction 

Hard scattering in QCD is accompanied by soft radiation of gluons caused by the drastic 
accelerations of color charge. The resulting antenna patterns are predicted to be particularly 
simple for 2^2 hard scattering processes in which one of the final particles is colorless § . 
The possibility of using the radiation pattern to identify the underlying parton structure of 
observed events has been called "partonometry" [Q]. In this paper, we address the question 
of how partonometry can be carried out in practice, and how it can be used to make further 
studies of QCD. 

As a specific example that is accessible to experiment, we study the process pp + jet 

at the Fermilab Tevatron energy y/s = 1.8 TeV. The predicted cross section for pj^* > 30 GeV/c 
with W ^ eu or fiu is ~ 180 pb, which is large enough for an ample number of events to be 
collected in the CDF and DO experiments at the Fermilab Tevatron. Similar physics can also 
be studied in or direct photon -|- jet production. 
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At lowest order in QCD, W^ + jet production comes from three incoherent partonic subpro- 
cesses: 

gq Wq (1) 
qg ^ Wq (2) 
qq^Wg (3) 

The first two of these are "Compton" subprocesses that are distinguished from each other by 
the direction of their initial gluon, which is in the +z direction for (|1|) and —z for (|^). The 
third "annihilation" subprocess has a slightly larger cross section than the sum of the Compton 
subprocesses, because of the dominance of quarks over gluons in the parton distributions at 
the moderate values of x involved here. 

"Partonometry" in this instance is the ability to recognize which of the three hard sub- 
processes (0)-(0) was responsible for a given W + jet event. In section 2 we will see that 
the subprocesses have very distinct radiation patterns. These patterns are quite striking, but 
unfortunately the ability to distinguish them is hampered by large fluctuations of the radiation 
around these average distributions. To model the fluctuations, as well as the underlying back- 
ground event and hadronization effects, we use the Monte Carlo program HERWIG 0, which 
has the correct correlations included to produce the radiation patterns that we seek. In section 
3 we define some event variables that are useful for distinguishing the underlying subprocesses. 
In section 4 we describe the Monte Carlo simulation that we use and show that it reproduces 
the radiation patterns that we expect, even after the underlying event and hadronization have 
been included. Our main results are in section 5, where we identify an observable of the 
soft radiation which can be used to discriminate between the parton-level subprocesses, and 
we show how this observable can be tested by experiment. Finally, in section 6 we give our 
conclusions. 



2. Radiation Patterns 

In the large Nc limit, soft gluon radiation from the three basic subprocesses is proportional to 

cosh /Xr] — cos A0 

cosh Arj — cos A0 

cosh Ar] 
cosh At] — cos A(f) 

for processes (1), (2), and (3), respectively. The "Lego" variables here are defined relative to 
the jet location: 

Art = 7]- r^jet (7) 
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Figure 1: Radiation densities (a) pi and (c) p^. p2 can be obtained from pi by the replacement 
At] — >■ —Arj. 

and 

A(j) = (j)- (^jet, (8) 

with 1] = — In tan ^/2 the pseudorapidity and (p the azimuthal angle. The soft gluon radiation 
patterns are displayed in cylindrical-coordinate plots in Fig. ^ Here At] is along the cylindrical 
axis, A0 is the azimuthal angle, and the intensity of the radiation is proportional to the distance 
from the cylindrical axis. Note that these radiation patterns only depend on the coordinates 
relative to the jet and are completely independent of the direction of the W boson. 

There are three important gross features of the radiation patterns. A first obvious feature 
is that the radiation patterns are singular in the direction coUinear to the outgoing jet. This 
is exhibited by the fact that the surfaces in Fig. |I| blow up to infinity as Arj — >■ 0, A(j) —>■ 0. 
In addition, the radiation from a gluon jet is stronger than radiation from g or g as expected: 
Pi= P2 = |P3 ^ 1/[(A^)^ + (A0)^] for small Ar], A0. 

The second general feature is that radiation from the Compton subprocesses is enhanced 
in the direction of their initial gluon: 

Pi ^ n Pi ^ 2 ] 

P2 ^ 2 Wor Ar; > 1 , ^ 1 > for Ar/ < 1 . (9) 

P3 1 J P3 ^ 1 J 

This enhancement extends arbitrarily far away from the jet in the [t], (p) plane. The large 
area over which the radiation pattern can be integrated offers the prospect of observing it in a 
majority of events, even though the observed radiation in any small region of [r], 0) is subject 
to large fluctuations and is partially masked by the "background event" . Note that we use the 
variables At] and A0 rather than AR and (3 of Ref. 0] because we have found regions distant 
from the jet to be important. The relations between these variables are Arj = AR cos (3 and 
A(p = AR sin (3. 

The third important gross feature of the radiation pattern is the suppression of radiation 
from the annihilation subprocess @ at small Ai] on the side opposite to the jet in azimuthal 
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Figure 2: (a) Soft gluon radiation function averaged over strip < \A(f)\ < 60°. Dashed curve 
= subprocess (|ip, Dotted curve = subprocess (0), Solid curve = subprocess (^. (b) Same, but 
for 120° < |A0| < 180° 



angle. The extreme of this is 

Pi 

P2 
P3 



1.25 1 

1.25 I for Ar] = 0, = 180' 
0.5 J 



(10) 



These features of the radiation patterns can be seen more quantitatively in Fig. 0, which 
shows the Arj distributions averaged over strips < |A0| < 60° and 120° < |A0| < 180°. 
Although we use the large Nc limit throughout, we note that the corrections are small 

and actually enhance the effects that we are seeking. 

3. Event Variables 

Before we look for soft radiation patterns in the data, we define some convenient observable 
quantities in order to highlight the features discussed above. We divide the region outside 
the jet into forward (F), backward (B) and side (S) regions as in Fig. |[ We then define the 
following scalar sums 

F= E b^h (11) 

+1.25 <Aj7< +3.75 

B= E bJ> (12) 

-3.75 <Ar?< -1.25 
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Figure 3: Definition of regions F, B and S. 

and 

s= E \P±\- (13) 

|A??|<1.25, \A(f>\>TT/2 

Witli tliese definitions, we would expect on average that subprocess (1) {gq — > Wq) would 
have an enhancement in F, subprocess (2) {qg — >■ Wq) would have an enhancement in B, and 
subprocess (3) {qq — >• Wg) would have a deficit in S. We have investigated a number of observ- 
able quantities that use this information to distinguish between the underlying subprocesses. 
The details are given in the appendix. In the main part of the text we will focus on the 
quantity \F — B\ + S . It is designed to be sensitive to the differences between events with a 
final-state quark (1 and 2) and events with a final-state gluon (3). This is because quark-jet 
events tend to have large forward-backward differences and relatively large radiation in the 
side region, as we saw in the previous section. Gluon-jet events, on the other hand, tend to 
have equal forward and backward radiation, in addition to suppression in the side region. 

To test this observable, it will be helpful to have another quantity sensitive to the dif- 
ferences between the various subprocesses that does not depend on the properties of the soft 
radiation. Various ways to measure the radiation inside a jet cone have been invented in order 
to discriminate between quark jets and gluon jets. They have been validated both by using 
Monte Carlo simulation [|, || and experiment 0. One suitable variable can be defined by 
dividing the jet cone of radius R = 0.7 in (r/, (p) into 0.1 x 0.1 cells, as is typical of calorimetric 
detectors. The jet p± can be defined as the total p± detected in all cells in the cone. By sorting 
these cells in order of their \p±\, we define the variable 

X = the minimum number of cells needed to include all but \J\p±\^°^ of the total 

(14) 

where p± is given in units of GeV/c. This definition makes the X distribution roughly inde- 
pendent of jet p±. X can be defined to include a fractional part given by linear or logarithmic 
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interpolation. 



4. Simulation 

The distinctive radiation patterns of Fig. |^ are predicted by the soft gluon approximation. 
However, on an event-by-event basis there will be large fluctuations around these average 
distributions. We must also ask to what extent these patterns survive finite-energy effects 
and the non-perturbative effects of hadronization, whereby colored partons combine to form 
physical hadrons. Many of those hadrons are resonances that undergo strong decay before 
reaching the detector, further broadening structures in the [r], 0) plane. The notion of local 
parton-hadron duality (LPHD) [^, which is supported by experiment suggests however 
that the patterns should survive reasonably well. 

For quantitative predictions, we simulate pp —>■ + jet using the HERWIG 5.8 |^ Monte 
Carlo event generator, which contains the color coherence effects we seek and embodies LPHD 
in its cluster formation and decay model. We make a rather tight cut on the jet transverse 
momentum (30 < j9j^* < 40 GeV/c) in order to simplify the study. We will also generally make 
a cut |?7jet| < 0.5 on the jet rapidity. These cuts are not essential, and it will be desirable to 
loosen them for actual data analysis to improve the statistics. In the experiment, the cut on 
minimum * should be replaced by a cut on pj^ to avoid biases that depend on jet structure. 
The effects we seek are somewhat cleaner at higher but the cross section falls rather 
rapidly at large 

Results from the simulation are shown in Fig. ^ for the same two strips in A0 as in Fig. |^. 
The quantity plotted is {p±) in GeV/c per unit Ar] Acf) at \i]jet\ < 0.5. A suppression at 
large t] is apparent, so that the "plateau" is not flat. But otherwise, the soft gluon pattern 
survives very well. In particular the roughly 2 : 1 enhancement of subprocesses (P and (0) 
in the directions of their incoming gluon at large At], and the suppression of subprocess @ 
at {rj ~ ?7jet, A0 ^ 180°} remain. In generating Fig. ^, the HERWIG option of a contribution 
from a minimum-bias-like "background event" was turned off. The result of turning it on is 
shown in Fig. ^. The background event adds extra radiation to each subprocess, so that the 
differences are largely preserved, while the ratios are substantially diluted. 

5. Observing the radiation patterns 

Each of the partonic subprocesses is predicted to have its distinctive radiation pattern, 

as we have seen. Unfortunately, these radiation patterns cannot be observed directly, because 
for real events — unlike the simulated ones — one does not know which of the three subpro- 
cesses is responsible for any given event. To establish this physics with a testable prediction, 
we must study a global correlation that is based on the fact that any particular event is caused 
by a single one of the three production subprocesses. Therefore, a typical event should have 
either enhanced radiation at ?7 ^ ?7jet; or enhanced radiation at r/ <^ rj^et] or reduced radiation 
at ?7 ~ r/jet, A(f) ^ 180°. Looking for these characteristics of the radiation in an event-by-event 
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Figure 4: (a) Average p± per unit area in for simulated events with |r7jet| < 0.5, in the 

region < |A0| < 60° (cf. Fig. D, "background event" off. (b) Same, but for 120° < |A0| < 
180°. 




V - Vjci 



(b) 



Figure 5: (a) Average p± per unit area in for simulated events with \r]jct\ < 0.5, in the 

region < |A</)| < 60° (cf. Fig. ||), "background event" on. (b) Same, but for 120° < |A0| < 
180°. 
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Figure 6: Probability distribution for \F — B\ + S. Solid curve: quark jets (from subprocess 
(|)); dotted curve: gluon jets (from subprocesses (|l|) and (^). 

basis is also what is needed to achieve the promise of "partonometry" Of course, in practice 
the underlying subprocess can only be determined with some probability, due to fluctuations 
in the soft radiation patterns. 

There are three necessary steps to partonometry. First, we must identify an observable 
function of the soft radiation that discriminates between partonic subprocesses. Second, we 
need to be able to test the prediction experimentally. Finally, we need a quantitative measure 
of the discrimination efficiency of our observable. The second step requires comparison with 
another event observable that does not depend on the properties of the soft radiation in the 
event. An observable based on information inside the jet cone is the natural choice. 

In what follows we give a concrete example of partonometry using the soft radiation observ- 
able \F — B\ + S . We also use the jet observable X as a cross check of this discriminator. These 
may or may not prove to be the best variables when detector effects are included. Nonetheless, 
the procedure should be applicable to a wide variety of observable quantities. 

The predictions from HERWIG for \F — B\ + S are shown in Fig. The differences between 
the curves satisfy our first criterion above. The predictions from HERWIG for the jet variable 
X are shown in Fig. ^ There is a clear difference between the quark jets (subprocess (H)) and 
the gluon jets (subprocesses (|1|) and (g)), with the latter tending toward larger values of X, 
i.e., gluon jets are fatter. The correlations between Figs. ^ and ^ provide the qualitative and 
quantitative tests for our soft radiation predictions. 

To observe the correlation, we divide the data sample into bins of \F — B\ + S and then 
calculate the mean value of X in each bin. The result is shown in Fig. ^. For a pure sample 
of gluon jets or quark jets, the value of (X) increases mildly with \F — B\ + S, indicating that 
(X) increases with the amount of soft radiation, as one would expect. However, we predict 
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Figure 7: Probability distribution for X defined as the minimum number of cells in the jet 
cone R < 0.7 required to contain all but ^J\p±\^°^ of the total jet in the cone, where 

ptot -g gjygj-^ GeV/c. Solid curve: quark jets (from subprocess (|])); dotted curve: gluon jets 
(from subpro cesses (P and (0)). 
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Figure 8: Mean value of X as a function of \F — B\ + S . 
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Figure 9: Fraction of jets which are gluon jets, function of \F — B\ + S. 

that the combined data will show a decrease in (X) with increasing \F — B\ + S. This is an 
unambiguous qualitative test of the soft radiation patterns. The decrease comes about because 
as \F — B\ + S increases, the likelihood that the outgoing jet is a quark increases, reducing 
the value of (X). This is shown in Fig. |, where we plot the fraction of outgoing gluon jets as 
a function of \F — B\ + S, as calculated using HERWIG. The gluon or quark jet fraction can 
be made as high as 80% by focusing on events with very large or small values of \F — B\ + S , 
albeit at the expense of statistics. Even stronger discrimination between the processes can be 
achieved by using both variables \F — B\ + S and X in tandem. 

6. Conclusions 

Soft radiation patterns in Ty + jet production provide a novel insight into the structure of 
QCD. In this paper we have identified the most important structures in the radiation patterns, 
and we have defined an observable, \F — B\ + S, that can be used to distinguish the three 
partonic subprocesses (|l|)-(0). We have also shown how this observable can be cross-checked 
against a second observable, (X), which relies only on information from inside the jet cone. 
The prediction, shown by the solid curve in Fig. ||, is a striking one: we predict a negative 
correlation between the amount of radiation inside the jet cone and the amount of radiation 
outside it. This arises solely because more radiation outside the cone is associated with a 
higher probability of the jet being a quark jet, as shown in Fig. |^. The negative correlation 
is otherwise opposite to the naive expectation that the radiation outside the cone should rise 
with the radiation inside - as is confirmed for quark jets alone or gluon jets alone in Fig. p. 
The ability to discriminate between the underlying soft radiation patterns is potentially useful 
as a tool for identifying and studying parton-level subprocesses. One could imagine using 
this information both as a means to separate the gluon from the quark parton distribution 



10 



functions in the initial state and as a method to study the differences between quark and gluon 
jets in the final state. 

We emphasize again that the procedure we have outlined for utilizing and testing the soft- 
radiation effects could be followed for a variety of soft-radiation and jet observables. Further 
work may be needed in particular to optimize the suggested observables to suit detector lim- 
itations. In addition, the details presented here rely on the use of HERWIG for producing the 
hadronization effects and the background event. Information on jet properties taken directly 
from experiment may be used to lessen this dependence on the Monte Carlo simulation. 
The important point is that the soft radiation effects should remain and can be discerned. Ob- 
serving these patterns is both interesting in its own right and as a novel probe of the underlying 
partonic structure of QCD. 
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Appendix 

In the text we have used the quantity \F — B\ + S to discriminate between the soft radiation 
patterns from the different event types. Many other possibilities are available. We give some 
of the alternatives here. 

One logical quantity to investigate is the asymmetry 

. F-B , , 



In the limit where the forward and backward regions extend to infinity, Eqs. predict the 

averages {A) = +1/3, — 1/3 and for subprocesses 1, 2, and 3, respectively. This prediction 
does not include such real- world effects as hadronization and the underlying event, which tends 
to increase both F and B. 

The forward-backward difference F — B itself is less sensitive to the underlying event. The 
histogram in Fig. p!0| shows the probability distribution for F — B. The peak is at F — i? > for 
subprocess (|I|), as expected from the average antenna patterns. The fraction of events of this 
type with F — B > is 67%. In other words, if one attempts to discriminate between equal 
numbers of subprocess (Q) and subprocess on the basis of F — B, the correct assignment is 
made 67% of the time. 

A second quantity that can be considered is just the total amount of soft radiation F+B+S. 
The histogram in Fig. |ll| shows the probability distributions for F + B + S for two different 
ranges of the rapidity of the outgoing jet. The distributions for 0.5 < jr^jetl < 2 are almost 
identical to those for more central jets \rijet\ < 0.5 , so in analyzing experiments, a large range 
of 77jet will be usable. 
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Figure 10: Normalized distributions oi F — B for events with |?7jet| < +0.5. Solid curve: quark 
jets (from subprocess (|I|)), dotted curve: gluon jets (from subprocess (0)). The distribution for 
quark jets from subprocess @ would be the same as for subprocess (|l]) with F — B replaced 
by-(F-5). 



Figure 11: Normalized distributions of F + + S". Dashed curve: quark jets (from subprocess 
(0) or (ID), dotted curve: gluon jets (from subprocess (§)). Normal dots and dashes are for 
|?7jet| < 0.5; sparse ones are for 1 < \rj-^et\ < 2. 
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The most important property of a soft-radiation observable is the abihty to distinguish 
events with outgoing gluon jets (subprocess @) from events with outgoing quark jets (subpro- 
cesses and (|^)). To quantify this abihty for a given observable O, we consider a sample of 
events that contains equal numbers of outgoing gluon and quark jets. This is close to the actual 
predicted mix. We then identify the 50% of events with the largest values of the observable O 
as quark-jet events and the 50% of events with the smallest values of O as gluon-jet events. 
The percentage of correctly-assigned events is then a measure of the discrimination power of 
the observable O at 100% efficiency. 

With this definition, we find a discrimination power for \F — B\ of 60%. As expected, 
the forward-backward difference is better at distinguishing subprocess (|l]) from subprocess 
than it is at distinguishing either from subprocess @. The total amount of radiation F + B + S 
does better, with a discrimination power of 63.6%. The quantity \F — B\ + S used in the body 
of the text improves upon this slightly further, having a discrimination power of 64.0%. We 
emphasize that these numbers are given at 100% efficiency. By rejecting events with middling 
values for the observable, one can increase the discrimination power at a cost of decreased 
efficiency. 

We have also explored other variations on the basic observables above, including: 

1. Ignoring contributions from cells of size 0.1 x 0.1 in (77, (p) that receive less than 0.1, 0.2, 
or 0.3 GeV/c of transverse momentum; 

2. Using multiplicity or charged multiplicity in place of the scalar sum of \p±\ as the measure 
of radiation intensity; 

3. Generalizing the observables to include n > 3 regions in (Ar^, A</)), each which may be 
assigned a different weight. 

These variations were found to be useless. In particular, variation (1) might have been expected 
to help by reducing the effect of the background event. On the contrary, however, provided 
HERWIG's model for the background event is correct, this is more than offset by a simultaneous 
reduction in the accuracy of measuring the signal. Variation (3) was explored extensively 
by using the Monte Carlo data itself. Cutting the (Ar/, Acp) plane into n = 3 — 8 discrete 
regions allowed each event to be represented by an ra-tuple vector consisting of the transverse 
momentum deposited into each of the n regions. One could then attempt to classify a given 
event, whose underlying subprocesses is taken to be unknown, by its closest resemblance to 
other events whose underlying subprocess is taken to be known, using a suitable metric. This 
procedure amounts to an idealization of what a neural network approach would only be able 
to approximate. However, it was not found to do any better than the classification based on 
F, B, and S. 

For conceptual simplicity, we have concentrated on the region —0.5 < r/jet < 0.5 . To 
increase statistics, this region can be extended to — 2 < r/jet < 2 with very little change in the 
symmetric distributions such as F + B + S or \F — B\+ S . Meanwhile, the ability to distinguish 
between subprocess (|I|) and subprocess on the basis of the antisymmetric quantity F — B 
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Figure 12: Distributions of jet rapidity: dashed curve = subprocess (|1]), dotted curve = 
subprocess (H), solid curve = subprocess (H). 

actually improves with increasing l^etl because of an accidental circumstance: the general 
radiation pattern peaks around 1] = for kinematical reasons. This tendency increasingly 
favors negative F — B as r^jet becomes positive; but subprocess (0) increasingly dominates over 
subprocess (|T]) in that limit because of parton distribution effects, as is shown in Fig. |T^. 
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